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Abstract. A Monte Carlo simulation has been performed to track light rays in 
cylindrical multimode fibres by ray optics. The trapping efficiencies for skew and 
meridional rays in active fibres and distributions of characteristic quantities for 
all trapped light rays have been calculated. The simulation provides new results 
for curved fibres, where the analytical expressions are too complex to be solved. 
The light losses due to sharp bending of fibres are presented as a function of the 
ratio of curvature to fibre radius and bending angle. It is shown that a radius of 
curvature to fibre radius ratio of greater than 65 results in a light loss of less than 
10% with the loss occurring in a transition region at bending angles <I> ~ 7r/8rad. 
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1. Introduction 

Active optical fibres are becoming more and more important in the field of detection 
and measurement of ionising radiation and particles. Light is generated inside the fibre 
either through interaction with the incident radiation (scintillating fibres) or through 
absorption of primary light (wavelength-shifting fibres). Plastic fibres with large core 
diameters, i.e. where the wavelength of the light being transmitted is much smaller 
than the fibre diameter, are commercially available and readily fabricated, have good 
timing properties and allow a multitude of different geometrical designs. The low costs 
of plastic materials make it possible for many present day or future experiments to 
use such fibres in large quantities (see ^ for a review of active fibres in high energy 
physics). For the construction of the highly segmented tracking detector of the ATLAS 
experiment approved for the LHC collider at CERN more than 600,000 wavelength- 
shifting fibres have been used |2]. Our work is also motivated by the fact that spiral 
fibres embedded in scintillators are being used for calorimetric measurements in long 
base line neutrino oscillation experiments, most recently in the MINOS experiment jSj- 

I Present address: Institut fiir Kernphysik, Joh. Gutenborg-Universitat Mainz, JJBecher-Weg 45, 
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The treatment of small diameter optical fibres involves electromagnetic theory 
applied to dielectric waveguides, which was first achieved by Snitzer 0] and Kapany 
ei aZ IS] • Although this approach provides insight into the phenomenon of total internal 
reflection and eventually leads to results for the field distributions and electromagnetic 
radiation for curved fibres, it is advantageous to use ray optics for applications to large 
diameter fibres where the waveguide analysis is an unnecessary complication. In ray 
optics a light ray may be categorised by its path along the fibre. The path of a 
meridional ray is confined to a single plane, all other modes of propagation are known 
as skew rays. The optics of meridional rays in fibres was developed in the 1950s [HI 
and can be found in numerous textbooks, e.g. in |7||Slin|- Since then, the scientific 
and technological progress in the field of fibre optics has been enormous. Despite the 
extensive coverage of theory and experiment in this field, only fragmentary studies on 
the trapping efficiencies and refraction of skew rays in curved multimode fibres could 
be found jlU[lll[[T^ll3| . We have therefore performed a three-dimensional simulation 
of photons propagating in simple circularly curved fibres in order to quantify the losses 
and to establish the dependence of these losses on the angle of the bend. We have also 
briefly investigated the time dispersion in fibres. For our calculations a common type 
of flbre in particle physics is assumed, specified by a polystyrene core of refractive 
index Ucore = 1-6 and a thin polymethylmethacrylate (PMMA) cladding of refractive 
index Udad — 1.49, where the indices are given at a wavelength of 590 nm. Another 
common cladding material is fiuorinated polymethacrylate with riciad = 1-42. Typical 
diameters are in the range of 0.5 - 1.5 mm. 

This paper is organised as follows: section 2 describes the analytical expressions of 
trapping efficiencies for skew and meridional rays in active, i.e. light generating, fibres. 
The analytical description of skew rays is too complex to be solved for sharply curved 
fibres and the necessity of a simulation becomes evident. In section 3 a simulation code 
is outlined that tracks light rays in cylindrical fibres governed by a set of geometrical 
rules derived from the laws of optics, section 4 presents the results of the simulations. 
These include distributions of the characteristic properties which describe light rays 
in straight and curved fibres, where special emphasis is placed on light losses due to 
the sharp bending of fibres. Light dispersion is briefly reviewed in the light of the 
results of the simulation. The last section provides a short summary. 

2. Trapping of photons 

When using scintillating or wavelength-shifting flbres in charged particle detectors the 
trapped light as a fraction of the intensity of the emitted light is very important in 
determining the light yield of the system. For very low light intensities as encountered 
in many particle detectors the photon representation is more appropriate to use than 
a description by light rays. Whether the flbres are scintillating or wavelength-shifting 
one is only ever concerned with a few lO's or lOO's of photons propagating in the fibre 
and single photon counting is often necessary. 

The geometrical path of any rays in optical fibres, including skew rays, was first 
analysed in a series of papers by Potter jj^ and Kapany |E1 ■ The treatment of angular 
dependencies in our paper is based on that. The angle 7 is defined as the angle of 
the projection of the light ray in a plane perpendicular to the axis of the fibre with 
respect to the normal at the point of refiection. One may describe 7 as a measure 
of the 'skewness' of a particular ray, since meridional rays have this angle equal to 
zero. The polar angle, 9' , is defined as the angle of the light ray in a plane containing 
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the fibre axis and the point of reflection with respect to the normal at the point of 
reflection. It can be shown that the angle of incidence at the walls of the cylinder, a, 
is given by cos a = cos 9' cos 7. The values of the two orthogonal angles 9' and 7 will 
be preserved independently for a particular photon at every reflection along its path. 

In general for any ray to be internally reflected within the cylinder of the fibre, 
the inequality sin a > aiii9'^^^^ = ridad/ncore must be fulfilled, where the critical angle, 
is given by the index of refraction of the fibre core, Ucore, and that of the cladding, 
ndad- In the meridional approximation the above equations lead to the well known 
critical angle condition for the polar angle, 9' > 9'^^,^^, which describes an acceptance 
cone of semi-angle, 9 [= 7t/2 — 9'], with respect to the fibre axis (see for example ^3] 
and references therein). Thus, in this approximation all light within the forward cone 
will be considered as trapped and undergo multiple total internal reflections to emerge 
at the end of the flbre. 

For the further discussion in this paper it is convenient to use the axial angle, 9, 
as given by the supplement of 9' , and the skew angle, 7, to characterise any light ray 
in terms of its orientation, see figure ^ for an illustration. 

The flux transmitted by a fibre is determined by an integration over the angular 
distribution of the light emitted within the acceptance domain, i.e. the phase space 
of possible propagation modes. Using the expression given by Potter et al ^HI and 
setting the transmission function, which parameterises the light attenuation, to unity 
the light flux can be written as follows: 

F — F A- F 

pScrit 

= V/ / / I{9,<P) cosi-fd^dn + 

Je=Q Ji=o J0=o UJ 

f.7r/2 pTT/2 pTT/2 

■2 1 I I Tin A\ „ 2 , 



4pM / / I{9, (/)) cos^ -fd-fdn , 

je=e„it Ji=~{e) J<p=o 

where dfl is the element of solid angle, j{9) refers to the maximum axial angle allowed 
by the critical angle condition, p is the radius of a cylindrical fibre and 1(9, (j)) is the 
angular distribution of the emitted light in the fibre core. The two terms, F„i and 
Fs, refer to either the meridional or skew cases, respectively. The lower limit of the 
integral for i^^ is 7 = arccos (sin ^crit/ sin^). 

The trapping efficiency for forward propagating photons, e^^^, may be defined 
as the fraction of totally internally reflected photons. The formal expression for the 
trapping efficiency, including skew rays, is derived by dividing the transmitted flux by 
the total fiux through the cross-section of the fibre core, Fq. For isotropic emission of 
fluorescence light the total flux equals Att^p^Iq. Then, the flrst term of equation (1) 
gives the trapping efficiency in the meridional approximation, 

e]/i = F™/Fo = i(l - cos 9,„t) « % , (2) 

where all photons are considered to be trapped if < 9 crit i independent of their actual 
skew angles. 

The integration of the second term of equation (1) gives the contributions of all 
skew rays to the trapping efficiency. Integrating by parts, one obtains 

^1/2 1 . „ cos^ 6* c„t sin 6* crit dt 



2 277 7o 7(1^(1 -tcos2^,„t) ' 
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with t = cos^ 9/ cos^ dcrit- Complex integration leads to the result: 

= i(l - cos 6'c„0 cosdcrit ■ 
The total initial trapping efficiency is then: 



(4) 



e^/'' ^ -{locos' 6 c„t) 




(5) 



which is approximately twice the trapping efficiency in the meridional approximation 
for small critical angles. The trapping efficiency of rays is crucially dependent on the 
circular symmetry of the core-cladding interface. Any ellipticity or variation in the 
fibre diameter will lead to the refraction of some skew rays. Furthermore, skew rays 
have a much longer optical path length, suffer from more reflections and therefore get 
attenuated more quickly, see section 4 for a quantitative analysis of this effect. In 
conclusion, for long fibres the effective trapping efficiency is closer to than to e^^^. 
Formula 121 yields a trapping efficiency of ell"^ — 3.44% for plastic fibres with Ucore = 
1.6 and Udad = 1.49. For these "standard" parameters the efficiency in formula 01 
evaluates to eV^ = 3.20% and in formulaElto e^/^ = 6.64%. 

3. Description of the tracking code 

The aim of the program is to track light rays through a fibred. Since the analytic 
analysis of the passage of skew rays along a curved fibre is exceedingly complex we 
treat the problem with a Monte Carlo technique. This type of numerical integration 
using random numbers is a standard method in the field of particle physics and is now 
practical given the CPU power currently available. On its path the ray is subject to 
attenuation, parameterised firstly by an effective absorption coefficient and secondly 
by a reflection coefficient. At the core-cladding interface the ray can be reffected 
totally or partially internally. In the latter case a random number is compared to the 
reflection probability to select reflected rays. 

Light rays are randomly generated on the circular cross-section of a fibre with 
radius p. An arbitrary ray is defined by its axial and azimuthal or skew angle. An 
advantage of this method is that any distribution of light rays can easily be generated. 
The axis of the fibre is defined by a curve z = f{s) where s is the arc length. For 
s < 0, it is a straight fibre along the negative 2:-axis and ioi < s < Lp, the fibre 
is curved in the x^-plane with a radius of curvature Rcurve- In particular, the curve 
/(s) is tangential to the z-axis at s = 0. 

Light rays are represented as lines and determined by two points, r and f* '. The 
points of incidence of rays with the core-cladding interface are determined by solving 
the appropriate systems of algebraic equations. In the case of a straight fibre the 
geometrical representation of a straight cylinder is used resulting in the quadratic 
equation 



The positive solution for the parameter m defines the point of incidence, r^, on the 
cylinder wall. In the case of a fibre curved in a circular path, the cylinder equation is 
generalized by the torus equation 



{x + (x' - a;) X rn)^ + {y + {y' - y) X- mf - = . 



(6) 



+ {z+{z' -z)xmf f'^f 

+{y + {y' ~y)y-mf - P^ = . 



(7) 
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The coefficients of this fourth degree polynomial are real and depend only on Rcurve 
and the vector components of r and r ' up to the fourth power. In most cases there 
are two real roots, one for the core-cladding intersection in the forward direction and 
one at r if the initial point already lies on the cylinder wall. The roots are found 
using Laguerre's method jl7| . It requires complex arithmetic, even while converging 
to real roots, and an estimate for the root to be found. The routine implements a 
stopping criterion in case of non-convergence because of round-off errors. The initial 
estimate is given by the intersection point of the light ray and a straight cylinder that 
has been rotated and translated to the previous reflection point. A driver routine is 
used to apply Laguerre's method to all four roots and to perform the deflation of the 
remaining polynomial. Finally the roots are sorted by their real part. The smallest 
positive, real solution for m is then used to determine the reflection point, r^j. 

After the point of incidence has been found, the reflection length and absorption 
probability can be calculated. The angle of incidence, a, is given by cos a = • n, 
where n denotes the unit vector normal to the core-cladding interface at the point of 
reflection and fin = {f — fn) / \r — fp\ is the unit incident propagation vector. Now the 
reflection probability corresponding to this angle a is determined. In case the ray is 
partially or totally internally reflected the total number of reflections is increased and 
the unit propagation vector after reflection, font, is then calculated by mirroring fm 
with respect to the normal vector: font ~ nn — 2n cos a. The program returns in a loop 
to the calculation of the next reflection point. When the ray is absorbed on its path or 
not reflected at the reflection point the next ray is generated at the fibre entrance end. 
A scheme on the main steps of the program can be found in figure |21 At any point of 
the ray's path its axial, azimuthal and skew angle are given by scalar products of the 
ray vector with the coordinate axes in a projection on a plane perpendicular to the 
fibre axis and parallel to the fibre axis, respectively. The transmitted flux of a specific 
fibre, taking all losses caused by bending, absorption and reflections into account, is 
calculated from the number of lost rays compared to the number of rays reaching the 
fibre exit end. 

This method gives rise to an efficient simulation technique for fibres with constant 
curvature. It is possible to extend the method for the study of arbitrarily curved 
fibres by using small segments of constant curvature. In the current version of the 
program light rays are tracked in the fibre core only and no tracking takes place in 
the surrounding cladding, corresponding to infinite cladding thickness. In long fibres 
cladding modes will eventually be lost, but for lengths < 1 m they can contribute 
to the transmission function. The simulation code is written in Fortran and it takes 
about 1.5 ms to track a skew ray through a curved fibre. 

4. Results of the tracking code 

Figuren]shows the passage of a skew ray along a straight fibre. The light ray has been 
generated off-axis with an axial angle oi 9 = 0.42 and would not be trapped if it were 
meridional. In general, the projection of a meridional ray on a plane perpendicular 
to the fibre axis is a straight line, whereas the projection of a skew ray changes its 
orientation with every reflection. In the special case of a cylindrical fibre all meridional 
rays pass through the fibre axis. The figure illustrates the preservation of the skew 
angle, 7, during the propagation. 



Computational studies of light in straight and curved fibres 6 
4.I. Trapping efficiency and acceptance domain 

Figure Ofa) shows the total acceptance domain and its splitting into the meridional 
and skew regions in the meridional ray approximation. The phase space density, i.e. 
the number of trapped rays per angular interval, is represented by proportional boxes. 
The density increases with cos^7 and sin 6. The contours relate to sharply curved 
fibres and are explained in section ^31 Figure l^tb) shows a projection of the phase 
space onto the sin^-axis. A peak around the value sinOcrit is apparent. A skew ray 
can be totally internally reflected at larger angles than meridional rays and the 
relationship between the minimum permitted skew angle, 7, at a given axial angle, 9, 
is determined by the critical angle condition: cos 7 = sin ^crit/ sin 0. 

Photons are generated randomly on the cross-section of the fibre with an isotropic 
angular distribution in the forward direction. The figure gives values for the two 
trapping efficiencies which can be determined by integrating over the two angular 
regions. The integrals are identical in value to the expressions in formulae |21 and 0] 
It is obvious from the critical angle condition that a photon emitted close to the 
cladding has a higher probability to be trapped than when emitted close to the centre 
of the fibre. For a given axial angle the range of possible azimuthal angles, in which 
the photons are uniformly distributed, for the photon to get trapped increases with 
the radial position, p, of the light emitter in the fibre core. It can be deduced 



from figure 4(a) that the meridional approximation is a good estimate for e if the 
photons originate at radial positions p < 0.8. The trapping of skew rays only becomes 
significant for photons originating at radial positions p > 0.9. This fact has been 



discussed before, e.g. in jTH]. Figure 4(b) shows the the trapping efficiency as a 



function of the axial angle. All photons with axial angles below Ocrit are trapped 
in the fibre, whereas photons with larger angles are trapped only if their skew angle 
exceeds the minimum permitted skew angle. It can be seen that the trapping efficiency 
falls off very steeply with the axial angle. 



4-2. Propagation of photons 

The analysis of trapped photons is based on the total photon path length per axial 
fibre length, P, the number of internal reflections per axial fibre length, 77, and the 
optical path length between successive internal reflections. In, where we follow the 
nomenclature of Potter and Kapany. It should be noted that these three variables are 
not independent as P = 77 x Z^. 

Figure El shows the distribution of the normalised path length, P{6), for photons 
reaching the exit end of straight and curved fibres of 0.6 mm radius. The figure also 
gives results for curved fibres of two different radii of curvature. The distribution of 
path lengths shorter than the path length for meridional photons propagating at the 
critical angle is almost flat. It can easily be shown that the normalised path length 
along a straight flbre is given by the secant of the axial angle and is independent of 
other flbre dimensions: P{0) ~ sec 6*. In case of the curved fibre the normalised path 
length of the trapped photons is less than the secant of the axial angle and photons 
on near meridional paths are refracted out of the fibre most. 

The distribution of the normalised number of reflections, r](9), for photons 
reaching the exit end of straight and curved fibres is shown in figure El Again, 
the figure gives results for curved fibres of two different radii of curvature. The 
number of reflections a photon experiences scales with the reciprocal of the flbre radius. 
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In the meridional approximation the normaUsed number of reflections is related by 
simple trigonometry to the axial angle and the fibre radius: r]m{0) = ta,nd/2p. The 
distribution of rj^m based on the distribution of axial angles for the trapped photons, 
is represented by the dashed line. The upper limit, riiOcrit), is indicated in the plot by 
a vertical line. The number of reflections made by a skew ray, 77^(6*), can be calculated 
for a given skew angle: ris{d) = rjm{0) / cos^ . It is clear that this number increases 
significantly if the skew angle increases. From the distributions it can be seen that in 
curved fibres the trapped photons experience fewer reflections on average. 

Figure Oa) shows the distribution of the reflection length, Ir{0), for photons 
reaching the exit end of fibres of radius p — 0.6 mm. The refiection length will scale 
with the fibre radius. The left figure shows Z_r(0) for four different over-all fibre 
lengths, Lp = 0.5, 1, 2 and 3 m, and the attenuation characteristics of the fibre is 
made apparent by the non-vanishing attenuation parameters used. Short reflection 
lengths correspond to long optical path lengths and large numbers of reflections. 
Because of the many reflections and the long total paths traversed, these photons 
will be attenuated faster than photons with larger reflection lengths. This reveals 
the high attenuation of rays with large skew angles. In the meridional approximation 
the refiection length is related to the axial angle by: l^ = 2p/ cosO. In the figure 
the minimum reflection length allowed by the critical angle condition is shown by 
a vertical line at lR{Ocrit) — 3.29mm. On average photons propagate with smaller 
reflection lengths along the curved fibre. Figure |7Jb) shows the distribution of Ir{6) 
in curved fibres of two different radii of curvature, Rcurve = 2 and 8 cm. The fibre 
radius is identical to the one used in figure (a) and the over-all fibre length is 0.5 m. 
It can be seen that the sharp peak in the photon distribution flattens with decreasing 
radius of curvature, so that the region of highest attenuation is close to the reflection 
length for photons propagating at the critical angle. 

In contrast to the analysis of straight fibres an approximation of the sharply 
curved fibre by meridional rays is not a very good one, since only a very small fraction 
of the light rays have paths lying in the bending plane. It is clear that when a fibre 
is curved the path length, the number of reflections and the reflection length of a 
particular ray in the fibre are affected, which is clearly seen in figures |S1 IHI and [7|[b) . 
The over-all fibre length for the curved fibres in these calculations is 0.5 m and the 
fibres are curved for their entire length. The average optical path length and the 
average number of reflections in a fibre curved over a circular arc are less than those 
for the same ray in a straight fibre for those photons which remain trapped. 

^.3. Light attenuation 

Light attenuation in active fibres has many sources, among them self-absorption, 
optical non-uniformities, refiection losses and absorption by impurities. The two main 
sources of attenuation in this type of fibres are the absorption of scintillation light and 
Rayleigh scattering from small density fiuctuations. The self-absorption is due to the 
overlap of the emission and absorption bands of the fluorescent dyes. The cumulative 
effect of these attenuation processes can be conveniently parameterised by an effective 
attenuation length over which the signal amplitude is attenuated to 1/e of its original 
value, a method often applied in high energy physics applications. The attenuation of 
active fibres at wavelengths close to its emission band (400 — 600 nm) is much higher 
than in wavelength regions of interest for standard applications of communication 
fibres where mainly infrared light is transmitted (0.8 — 0.9 /im and 1.2 — 1.5 /im). 



Computational studies of light in straight and curved fibres 



8 



Restricting the analysis to these processes, the transmission through an active 
fibre can be represented for any given axial angle by T = ex-p[—P{9)Lp/Ai,uik] x 
qri{9)LF ^ where the exponential function describes light losses due to bulk absorption 
and scattering (bulk absorption length Af,uik)i and the second factor describes light 
losses due to imperfect reflections (reflection coefficient q) which can be caused by 
a rough surface or variations in the refractive indices. A comparison of some of our 
own measurements to determine the attenuation length of plastic fibres with other 
available data indicates that a reasonable value for the bulk absorption length is 
Abuik ~ 3 m. Most published data suggest a deviation of the reflection coefficient, 
which parameterises the internal reflectivity, from unity between 5 x 10~^ and 
6.5 X 10~^ ^ni- A reasonable value oi q = 0.9999 is used in the simulation to account 
for all losses proportional to the number of reflections. 

Internal reflections being less than total give rise to so-called leaky or non-guided 
modes, where part of the electromagnetic energy is radiated away. Rays in these modes 
populate a region defined by axial angles above the critical angle and skew angles 
slightly larger than the ones for totally internally reflected photons. These modes are 
taken into account by using the Fresnel equation for the refiection coefficient, (R), 
averaged over the parallel and orthogonal plane of polarisation 

/m_l^p IP ^ 1 / tan2(a-/3) , sin^(a~/3) \ 

where a is the angle of incidence and /3 is the refraction angle. However, it is obvious 
that non-guided modes are lost quickly in a small fibre. This is best seen in the 
fraction of non-guided to guided modes, /, which decreases from / = 11% at the first 
reflection of the ray over / = 2.5% at the second reflection to / < 1% at further 
reflections. Since the average reflection length of non-guided modes is Ir ~ 1.5 mm 
those modes do not contribute to the flux transmitted by fibres longer than a few 
centimeters. The absorption and emission processes in fibres are spread out over a 
wide band of wavelengths and the attenuation is known to be wavelength dependent. 
For simplicity only monochromatic light is assumed in the simulation and highly 
wavelength-dependent effects like Rayleigh scattering are not included explicitly. 

A question of practical importance for the estimation of the light output 
of a particular fibre application is its transmission function. In the meridional 
approximation and substituting exp(— Ing) by exp(l — q) the attenuation length can 
be written as 

A,n ^ cos 6 crit [i/Mulk + {'^ - q) Sine crtt/2p]~^ . (9) 

Only for small diameter fibres {D ~ 0.1mm) are the attenuation lengths due to 
imperfect reflections of the same order as the absorption lengths. Because of the large 
radii of the fibres discussed refiection losses are not relevant for the transmission 
function and the attenuation length contracts to A„j — A tuik cos 9 crit- For the 
simulated bulk absorption length this evaluates to A„i = 2.8 m. The transmission 
function outside the meridional approximation can be found by integrating over the 
normalised path length distribution, where dN represents the number of photons per 
path length interval dP, weighted by the exponential bulk absorption factor: 

1 f°° 

T=— dN/dPe-^^"'^"-"' dP . (10) 
^ Jp=o 

FigurelHlshows this transmission function versus the ratio of fibre to absorption length, 
Lp/ Am- A simple exponential fit, T oc exp [—Lp/Aejj], applied to the simulated light 
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transmissions for a varying fibre length results in an effective attenuation length of 
Aejff = 2.4 m. For Lp/Am > 0.2 this description is sufficiently accurate to parameterise 
the transmission function, at smaller values for Lp/Am the light is attenuated faster. 
The difference of order 15% to the meridional attenuation length is attributed to the 
tail of the path length distribution. 

Measurements of the light attenuation in fibres proves this simple model of a single 
attenuation length to be wrong. A dependence of the attenuation length with distance 
usually is observed [201 ■ The most important cause of this effect is the fact that the 
short wavelength components of the scintillation light is dominantly absorbed. This 
leads to a shift of the average wavelength in the emission spectrum towards longer 
wavelengths and to an increase in the effective attenuation length. 

4- 4- Trapping efficiency and transition losses in sharply curved fibres 

One of the most important practical issues in implementing optical fibres into compact 
particle detector systems are macro-bending losses. In general, some design parameters 
of fibre applications, especially if the over-all size of the detector system is important, 
depend crucially on the minimum permissible radius of curvature. By using waveguide 
analysis transition and bending losses have been thoroughly investigated and a loss 
formula in terms of the Poynting vector can be derived |21l 122] ■ Those studies are 
difficult to extend to multimode fibres since a large number of guided modes has to 
be considered. When applying ray optics to curved multimode fibres the use of a 
two-dimensional model is common 

[nmSlini- in contrast, our simulation method 
follows a three-dimensional approach. 

Photons are lost from a fibre core both by refraction and tunnelling. In the 
simulation only refracting photons were considered. The angle of incidence of a light 
ray at the tensile (outer) side of the fibre is always smaller than at the compressed side 
and photons propagate either by reflections on both sides or in the extreme meridional 
case by reflections on the tensile side only. If the fibre is curved over an arc of constant 
radius of curvature photons can be refracted, and will then no longer be trapped, at 
the very first reflection point on tensile side. Therefore, the trapping efficiency for 
photons entering a curved section of fibre towards the tensile side is reduced most. 
Figure |51 quantifies the dependence of the trapping efficiency on the azimuthal angle, 
between the bending plane and the photon path for a curved fibre with a radius 
of curvature Rcurve = 2 cm. The azimuthal angle is orthogonal to the axial angle and 
XE* = rad corresponds to photons emitted towards the tensile side of the fibre. 

Figure [TUl displavs the explicit dependence of the transmission function for fibres 
curved over circular arcs of 90° on the radius of curvature to fibre radius ratio for 
different fibre radii, p = 0.2, 0.6, 1.0 and 1.2 mm. No further light attenuation is 
assumed. Evidently, the number of photons which are refracted out of a sharply 
curved fibre increases very rapidly with decreasing radius of curvature. The losses are 
dependent only on the curvature to fibre radius ratio, since no inherent length scale 
is involved, justifying the introduction of this scaling variable. The light loss due to 
bending of the fibre is about 10% for a radius of curvature of 65 times the fibre radius. 

The use of the meridional approximation in the bending plane in place of a three 
dimensional fibre is justified by the losses being dominantly caused by meridional 
rays ■ FigureEJshows a section of a curved fibre and the passage of a meridional 

ray in the bending plane with maximum axial angle. In this model photons are guided 
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for axial angles 

cos 6* < cosOo = — — — -cosOcrit , (11) 
R + p 

where the subscript O refers to the outer fibre wall. A transmission function can be 
estimated by assuming that all photons with axial angles > Oq are refracted out of 
the fibre: 

„ 1 cos Orrii 

T =1 . (12) 

1 + Rcurve/ P 1 ^ COS 9 crit 

This transmission function is shown in figure [Till as a dashed line and it overestimates 
the light losses due the larger axial angles allowed for skew rays. A comparable 
theoretical calculation using a two-dimensional slab model and a generalized Fresnel 
transmission coefficient has been performed by Badar et al ^2]- Their plot of the 
power contained in the fibre core as a function of the radius of curvature (figure 5) is 
similar to our results on the transmission function in the meridional approximation. 
In |10j a ray optics calculation for curved multimode fibres involving skew rays is 
presented. In this paper a discussion on the transmission function is missing. Instead, 
a plot of the power remaining in a curved fibre versus distance is shown which gives 
complementary information. 

For photons entering a curved section of fibre the first point of reflection on the 
tensile side defines the transition angle, trans i measured from the plane of entry. The 
angular range of transition angles associated with each ray is called the transition 
region of the fibre. For photons emitted towards the tensile side the transition angle is 
related to the axial angle and since the angular phase space density of trapped photons 
is highest close to the critical angle a good estimate is trans = 9 crit — do- For a fibre 
radius p = 0.6 mm and radii of curvature Rcurve = 1, 2, and 5 cm the above formula 
leads to transition angles '^trans — 0.19, 0.08 and 0.03 rad, respectively. We attribute 
these discrete angles to beams emitted in the bending plane. Photons emitted from 
the fibre axis towards the compressed side are not lost at this side, however, they 
experience at least one reflection on the tensile side if the bending angle exceeds the 
limit $iimit = a,T:ccos[Rcurve I {Rcurve + 2/9)] K. arccos [1 - 2 pI Rcurve]- A transition in 
the transmission function should occur at bending angles between ^umit/^, where all 
photons emitted towards the tensile side have experienced a reflection, and ^Hmit, 
where this is true for all photons. Figure [T^ shows the transmission as a function of 
bending angle, $, for a standard fibre as defined before. Once a sharply curved fibre 
with a ratio Rcurve/ P > 83 is bent through angles $ ~ 7r/8rad hght losses do not 
increase any further. The transition region ranges from 0.44 to ^ 1.06 rad and is 
indicated in the figure by arrows. At much smaller ratios Rcurve/ P the model is no 
longer valid to describe this behaviour. 

Experimental results on losses in curved multimode fibres along with 
corresponding predictions are best known for silica fibres with core radii p « 50 pm. 
Calculations on the basis of ray optics for a plastic fibre with p — 0.49 mm can be 
found in (13( . Our result on the transmission function in the meridional approximation 
T = 0.35 at p/ Rcurve =20 is in good agreement with the two-dimensional calculation. 
The larger value of T = 0.65 predicted by the simulation is explained by the small loss 
of skew rays, clearly seen in figured It should be noted that the difference between 
finite and infinite cladding and the appearance of oscillatory losses in the transition 
region has not been investigated in the simulation. 
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Figure 01 shows contours of the angular phase space for photons which were 
trapped in the straight fibre section but are refracted out of sharply curved fibres 
with radii of curvature Rcurve = 2 and 5 cm. The contours demonstrate that only 
skew rays from a small region close to the boundary curve are getting lost. The 
smaller the radius of curvature, the larger the affected phase space region. 

4-5. Light dispersion 

The timing resolution of scintillators are often of paramount importance, but a pulse 
of light, consisting of several photons propagating along a fibre, broadens in time. In 
active fibres, three effects are responsible for the time distribution of photons reaching 
the fibre exit end. Firstly the decay time of the fluorescent dopants, usually of the 
order of a few nanoseconds, secondly the chromatic dispersion in a dispersive medium, 
and thirdly the fact that photons on different paths have different transit times to reach 
the fibre exit end, known as inter-modal dispersion. 

The chromatic dispersion is due to the spectral width, AA, of the emitter. It is 
the combination of material dispersion and waveguide dispersion. If the core refractive 
index is explicitly dependent on the wavelength, n(A), photons of different wavelengths 
have different propagation velocities along the same path, called material dispersion. 
The broadening of a pulse is given by At = Lp/ccore {^^(Pn/dX^) AX/X 0. 
The FWHM of the emission peaks of scintillating or wavelength-shifting fibres is 
approximately 40 — 50 nm. The material dispersion in the used polymers (mostly 
polystyrene) is of the order of ns/nmxkm and thus negligible for multimode fibres. 

The transit time in ray optics is simply given by r = P{9) / Ccore, where Ccore is the 
speed of light in the fibre core. The simulation results on the transit time are shown 
in figure El The full widths at half maximum (FWHM) of the pulses in the time 
spectrum are presented for four different fibre lengths. The resulting dispersion has to 
be compared with the time dispersion in the meridional approximation which is simply 
the difference between the shortest transit time t{0 — 0) and the longest transit time 
t{9 = Ocrit)'- At = Lp/ Ccore {sec 6 crit — 1), whcrc Lp is the total axial length of the 
fibre. The dispersion evaluates for the different fibre lengths to 197 ps for 0.5 m, 393 ps 
for 1 m, 787 ps for 2 m and 1181 ps for 3 m. Those numbers are in good agreement with 
the simulation, although there are tails associated to the propagation of skew rays. 
With the attenuation parameters of our simulation the fraction of photons arriving 
later than t{9 = 9crit) decreases from 37.9% for a 0.5m fibre to 32% for a 3m fibre due 
to the stronger attenuation of the skew rays in the tail. Due to inter-modal dispersion 
the pulse broadening is quite significant. 

5. Summary 

We have simulated the propagation of photons in straight and curved optical fibres. 
The simulations have been used to evaluate the loss of photons propagating in fibres 
curved in a circular path in one plane. The results show that loss of photons due to 
the curvature of the fibre is a simple function of radius of curvature to fibre radius 
ratio and is < 10% if the ratio is > 65. The simulations also show that for larger ratios 
this loss takes place in a transition region ($ ^ during which a new distribution 
of photon angles is established. Photons which survive the transition region then 
propagate without further losses. 
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We have also used the simulation to investigate the dispersion of transit times of 
photons propagating in straight fibres. For fibre lengths between 0.5 and 3 m we find 
that approximately two thirds of the photons arrive within the spread of transit times 
which would be expected from the use of the simple meridional ray approximation 
and the refractive index of the fibre core. The remainder of the photons arrive in a 
tail at later times due to their helical paths in the fibre. The fraction of photons in 
the tail of the distribution decreases only slowly with increasing fibre length and will 
depend on the attenuation parameters of the fibre. 

We find that when realistic bulk absorption and reflection losses are included in 
the simulation for a straight fibre, the overall transmission can not be described by 
a simple exponential function of propagation distance because of the large spread in 
optical path lengths between the most meridional and most skew rays. 

We anticipate that these results on the magnitude of transition losses will be of 
use for the design of particle detectors incorporating sharply curved active fibres. 
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(a) (b) 

Figure 1. (a) Definition of tiic angles for a skew ray from an initial point, r to a 
point of reflection, rjj, on the core-cladding interface. The axial angle of the ray 
is 6. The angle between the projection of the ray in a plane perpendicular to the 
axis of the fibre with respect to the normal at the point of reflection, 7, defines 
the ray's 'skewness'. The angle of incidence at the point of reflection is a. (b) 
The helical path of a skew ray in a fibre of radius p = 0.6 mm which spirals along 
the core-cladding interface. The z-axis represents the axial fibre length and the 
fibre cross-section lies in the xy-plane. 
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Define fibre pargimeters: 
- Stra ight and bent section length, p, 
LOOP over photons 



Generate one photon: 

- Posi tion cind angle distribution according to emitter type 
OOP over reflections at core-cladding interface 
Calculate photon parameters at point r: 
- axial and azimuthal angle, skew ein gle 



Find next reflection point ra: 
Phot on absorbed on path? 

(yes 

Phot on reached fibre end face? 

^ YES 

Photon reflected at fji? 

Calculate propagation parameters: 

- reflection length, total path length, no 

Coordinate transformation fii—>-f 



numerical methods 



of reflections 



RETURN 



RETURN 



Calculate flux parameters: 

- bending, absorption eind reflection losses, trapping efficiency 



Figure 2. A scheme on the main steps of the program to generate and track 
photons in fibres of constant curvature. The equations to find the intersection of 
the Ught ray with the core-cladding interface are solved by numerical methods. 
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Figure 3. (a) The angular phase space domain for trapped photons in a fibre. 
To the left of the dividing line at sin dcrit all skew angles are accepted. To the 
right of the line a minimum skew angle is required by the critical angle condition. 
The trapping efficiencies are evaluated by integration over the two regions. The 
label 'Losses' points to contours for photons refracted out of sharply curved fibres 
with radii of curvature Rcurvc = 2 and 5 cm. (b) A projection of the phase space 
onto the sinS-axis. 
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Figure 4. Trapping efficiency, e^^^, for pliotons propagating in the forward 
direction as a function of radial position, p, of the hght emitter in the fibre core 
(a) and of the axial angle (b). The maximum axial angle allowed by the critical 
angle condition in the meridional approximation is indicated. 
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Figure 5. The distribution of the optical path length, P{9), of trapped photons 
in fibres of radius p = 0.6 mm normalised to the axial length of the fibre. The 
figure shows P{0) for a straight fibre and for two different radii of curvature, 
Rcurve = 2 and 8cm. The vertical line at P{Ocrit) = 1-074 indicates the upper 
limit of P in the meridional approximation. 
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Figure 6. The distribution of the number of reflections, r]{6), for trapped photons 
in fibres of radius p = 0.6 mm normalised to the axial length of the fibre. The 
figure shows ri{6) for a straight fibre and for two different radii of curvature, 
Rcurve = 2 and 8cm. The vertical line at r)m(dcrit) = 0.326 indicates its upper 
limit in the meridional approximation. The dashed line shows the distribution of 
rim{e) = tan0/2p. 




Figure 7. The distribution of the reflection length, Ir{6), for photons reaching 
the exit end of fibres of radius p = 0.6 mm. The figure shows lii{0) for straight 
fibres (a) with three different fibre lengths, Lp = 0.5, 1, 2 and 3 m and for curved 
fibres (b) of 0.5 m length with radii of curvature Rcurve = 2 and 8 cm including 
a comparison to the straight fibre. In both figures the vertical line at lii{9crit) = 
3.29 mm indicates the lower limit of Ir in the meridional approximation. 
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Ratio of Fibre to Absorption Length 

Figure 8. Transmission function for a straiglit fibre witli a bulk light absorption 
length Aj,„(;t = 3 m and a reflection coeflicient q = 0.9999. The transmission as a 
function of the ratio of fibre to absorption length, Lp/Am, is calculated from the 
optical path length distribution. A simple exponential fit results in an effective 
attenuation length of A^,^ = 2.4 m. The dashed line shows the transmission 
function in the meridional approximation with A 
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Figure 9. Trapping ofliciency, e^/^, for photons in a sharply curved fibre with 
radius of curvature Rcurve = 2 cm and fibre radius p = 0.6 mm. The trapping 
efficiency is plotted as a function of the azimuthal angle, between the photon 
path and the bending plane, so that ^ = Orad corresponds to photons emitted 
towards the tensile side of the fibre. 
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Figure 10. The transmission function for fibres curved over a circular arc of 
90 ° is plotted as a function of the radius of the curvature to fibre radius ratio 
for different fibre radii, p = 0.2, 0.6, 1.0 and 1.4 mm. The dashed line is a simple 
estimate from the meridional approximation. 




Figure 11. Section of a curved fibre with radius p and radius of curvature 
Rcurve- The passage of a meridional ray in the bending plane with maximum 
axial angle is shown. 
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Figure 12. Transmission function for a curved fibre of radius p = 0.6 mm witii 
three different radii of curvature, Rcurve = 1, 2 and 5 cm, corresponding to tfie 
ratios Rcurve/ P = 17, 33 and 83, respectively. The ordinate is the fraction of 
photons reaching the fibre exit end as a function of the bending angle, ^, and the 
arrows indicate the transition region in the meridional approximation. 




Figure 13. The distribution of the transit time in nanoseconds for photons 
reaching the fibre exit end. For the fibre lengths Lp = 0.5, 1, 2 and 3 m the pulse 
dispersion (FWHM) of the transit time distribution is 200, 400, 775, and 1200 ps, 
respectively. 



